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1 Introduction 


Supersymmetry (SUSY) is a landmark in our efforts to construct a unified theory of all 
fundamental interactions observed in nature. At very high energies, close to the Planck 
scale (Mp) it is indispensable in constructing consistent string theories, and at low energies 
(~ 1 TeV) it seems unavoidable if the gauge hierarchy problem is to be resolved. Such a 
resolution provides a measure of the supersymmetry breaking scale Msusy ~ 0(1 TeV). 
There is indirect evidence for such a low-energy supersymmetry breaking scale, from the 
unihcation of the gauge couplings 0] and from lightness of the Higgs boson as determined 
from precise electroweak measurements, mainly at LEP |^. Furthermore, such a low 
energy SUSY breaking scale is also favored cosmologically. As is well known, i?-parity 
conserving SUSY models, contain in the sparticle spectrum a stable, neutral particle, 
identihable with the lightest neutralino (y), referred to as the LSP P]. It is important 
that such a LSP with mass, as low-energy SUSY entails, in the 100 GeV — 1 TeV region, 
may indeed provide the right form and amount of the highly desirable astrophysically and 
cosmologically Dark Matter (DM). The latest data about Cosmic Microwave Background 
(CMB) radiation anisotropies not only favour a flat (A; = 0 or Dq = 1), inflationary 
Universe, but they also determine a matter density » 0.15 ± 0.05. Taking into 

account the simultaneously determined baryon density Dphg ^ 0.02, and the rather tiny 
neutrino density, they result to 


= 0.13 ±0.05. (1) 

If we assume that all DM is supersymmetric due to LSP, i.e. ^Idm = h is tempt¬ 
ing to combine the bound of Eq.^ with other presently available constraints from particle 
physics, such as the lower bound on the mass of the Higgs bosons {rrih > 113.5 GeV) 
provided by LEP and the recent results from the BNL E821 experiment on the 
anomalous magnetic moment of the muon (^a^ = 43(16) x 10“^°). Although the situ¬ 
ation regarding the — 2 has not been dehnitely settled, supersymmetry emerges as a 
prominent candidate in explaining the discrepancy between the Standard Model predic¬ 
tions and experimental measurements, and in the sequel we concede that this deviation 
accounted for SUSY. We hnd that this combination of the experimental information from 
high energy physics and cosmology puts austere bounds on the parameter space of the 
Constrained Minimal Supersmmetric Standard Model (CMSSM), enabling us to investi¬ 
gate the potential of discovering SUSY, if it is based on CMSSM, at future colliders and 
direct DM search experiments. 


2 Neutralino relic density 

It has been argued that for large tan (3 the neutralino relic density (D^ hg) can be com¬ 
patible with the recent cosmological data which favour small values for hg. In this 
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regime the neutralino (x) pair annihilation throngh s-channel psendo-scalar Higgs boson 
(H) exchange leads to an enhanced annihilation cross sections reducing signihcantly the 
relic density [0, while the heavy CP-even Higgs {H) exchange is P-wave suppressed and 
not that important. The importance of this mechanism, in conjunction with the recent 
cosmological data which favour small values of the DM relic density, has been stressed 
in [^,^. The same mechanism has been also invoked where it has been shown that it 
enlarges the cosmologically allowed regions. In fact cosmology does not put severe upper 
bounds on sparticle masses, and soft masses can be in the TeV region, pushing up the 
sparticle mass spectrum to regions that might escape detection in future planned accel¬ 
erators. Such upper bounds are imposed, however, by the recent ( 7 ^ — 2 E821 data 


constraining the CMSSM in such a way that supersymmetry will be accessible to LHC 
or other planned linear colliders if their center of mass energy is larger than about 
1.2 TeV 1^. The bounds put by the — 2 has been the subject of intense phenomeno¬ 
logical study the last few months |1T2|-[T7| . 


The XX fusion to the pseudo-scalar Higgs boson. A, which subsequently decays to a bb 
or a rr, becomes the dominant annihilation mechanism for large tan /?, when the pseudo- 
scalar mass rriA approaches twice the neutralino mass, rriA — 2m^. In fact by increasing 
tan/? the mass rriA decreases, while the neutralino mass remains almost constant, if the 
other parameters are kept hxed. Thus ttia is expected eventually to enter into the regime 
in which it is close to the pole value = 2m^, and the pseudo-scalar Higgs exchange 
dominates. It is interesting to point out that in a previous analysis of the direct DM 
searches 0, we had stressed that the contribution of the CP-even Higgs bosons exchange 
to the LSP-nucIeon scattering cross sections increases with tan /?. Therefore in the large 
tan (3 regime one obtains the highest possible rates for the direct DM searches. Similar 
results are presented in Ref. |^. In the framework of the CMSSM the chargino mass 
bound as well as the recent LEP Higgs mass bound ^ already exclude regions in which x 
has a large Higgsino component, and thus in the regions of interest the x is mainly a bino. 
A bino is characterized by a very small coupling to the pseudo-scalar Higgs A, however 
the largeness of tan (3 balances the smallness of its coupling giving a sizeable effect when 
rriA — 2m^, making the s-channel pseudo-scalar exchange mechanism important. 

It becomes obvious from the previous discussion that an unambiguous and reliable 
determination of the A-mass, m^, is necessary in order to to calculate the neutralino relic 
density especially in the large tan/3 region. The details of the procedure in calculating 
the spectrum of the CMSSM can be found elsewhere [^,|^. Here we shall only briefly 
discuss some points which turn out to be essential for a correct determination of m^. In 
the constrained SUSY models, such as the CMSSM, uia is not a free parameter but is 
determined once mo, M 1 / 2 , A as well as tan/? and the sign of fi are given. rriA depends 
sensitively on the Higgs mixing parameter, m|, which is determined from minimizing 
the one-loop corrected effective potential. For large tan /? the derivatives of the effec¬ 
tive potential with respect the Higgs helds, which enter into the minimization conditions, 
are plagued by terms which are large and hence potentially dangerous, making the per- 
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turbative treatment untrustworthy. In order to minimize the large tan/3 corrections we 
had better calculate the effective potential using as reference scale the average stop scale 
~ EH- this scale these terms are small and hence perturbatively valid. 


Also for the calculation of the pseudo-scalar Higgs boson mass all the one-loop correc¬ 
tions must be taken into account. In particular, the inclusion of those of the neutralinos 
and charginos yields a result for that is scale independent and approximates the pole 
mass to better than 2% 0. A more significant correction, which drastically affects the 
pseudo-scalar mass arises from the gluino-sbottom and chargino-stop corrections to the 
bottom quark Yukawa coupling hb The proper resummation of these corrections 

is important for a correct determination of hb |^,@, and accordingly of the mu- 


In calculating the LSP relic abundance, we solve the Boltzmann equation numerically 
using the machinery outlined in Ref. |^. In this calculation the coannihilation effects, in 
regions where tr approaches in mass the LSP, which is a high purity Bino, are properly 
taken into account. 

In what follows only the /i > 0 case is considered. The /r < 0 case is not favored by 
the recent 6 —>• sy data, as well as by the observed discrepancy of the — 2, if the latter 
is attributed to supersymmetry, and therefore we shall discard it. 

In the panels shown in figure |I] we display our results by drawing the cosmologically 
allowed region 0.08 < 0.18 (dark green) in the mo, M 1/2 plane for values of tan/3 

equal to 50 and 55 respectively. Also drawn (light green) is the region 0.18 <Vt^hQ< 0.30. 
The default values for the masses of massive quarks are rrit = 175 GeV, m,- = 1.777 GeV 
and mb{mb) = 4.25 GeV. The remaining inputs are shown on the top of each panel. The 
solid red mark the region within which the supersymmetric contribution to the anomalous 
magnetic moment of the muon falls within the E821 range 01 ^^^ = (43.0 ± 16.0) x 10“^°. 
The dashed red line marks the boundary of the region when the more relaxed lower 
bound 11 X 10“^*’ < is used, corresponding to the 2cr lower bound of the E821 

range. Along the blue dashed-dotted lines the light GP-even Higgs mass takes values 
113.5 GeV (left) and 117.0 GeV (right) respectively. The line on the left marks therefore 
the recent LEP bound on the Higgs mass 0. Also shown is the chargino mass bound 
104 GeV [|. The shaded area (in red) at the bottom of each figure, labelled by TH, is 
theoretically disallowed since the light stau is lighter than the lightest of the neutralinos. 
From the displayed figures we observe that for values of tan /3 up to 50 the cosmological 
data put an upper bound on the parameter uiq. However, there is practically no such 
upper bound for the parameter M 1 / 2 , due to the coannihilation effects [|T^] which allow 
for Mi /2 as large as 1700 GeV within the narrow coannihilation band lying above the 
theoretically disallowed region. 

For tan /3 = 55 a large region opens up within which the relic density is cosmologically 
allowed. This is due to the pair annihilation of the neutralinos through the pseudo- 


Tn the context of our analysis focus point regions show up for smaller values of the top mass. In 
any case the bulk of the focus point region appears for rather large values of mo and hence they are not 
favoured by the g^, — 2 data. 
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Figure 1: Cosmologically allowed regions of the relic density for two different values of 
tan/3 in the (Mi/ 2 ,rno) plane. The remaining inputs are shown in each figure. The mass 
of the top is taken 175 GeV. In the dark green shaded area 0.08 <Vt^h‘^< 0.18. In the 
light green shaded area 0.18 < /iq < 0.30 . The solid red lines mark the region within 
which the supersymmetric contribution to the anomalous magnetic moment of the muon 
is = (43.0 ± 16.0) x 10“^*^. The dashed red line is the boundary of the region for 

which the lower bound is moved to 11 x 10 “^° < The dashed-dotted blue lines 

are the boundaries of the region 113.5 GeV < mniggs < 117.0 GeV. 


scalar Higgs exchange in the s-channel. As explained before, for such high tan (3 the ratio 
mA/‘2m^ approaches unity and the pseudo-scalar exchange dominates yielding large cross 
sections and hence small neutralino relic densities. In this case the lower bound put by 
the Qg — 2 data cuts the cosmologically allowed region which would otherwise allow for 
very large values of mo, M 1 / 2 . The importance of these corridors has been stressed in the 
analysis of [^. However, in the analysis presented here these show up at higher values of 


the parameter tan (3. We should remark at this point that in our analysis we use the value 
of astrong{Mz) as iuput and relax unihcation of the 03 gauge coupling with the others. 
In the constrained scenario it is almost impossible to reconcile gauge coupling unihcation 
with a value for ctstrongiMz) consistent with experiment due to the low energy threshold 
effects. This change affects drastically the values of other parameters and especially that 
of the Higgsino (/i) and Higgs (m|) mixing parameters that in turn affect the pseudo-scalar 
Higgs boson mass which plays a dominant role. 

For the tan/3 = 55 case, close the highest possible value, and considering the con¬ 
servative lower bound on the muon’s anomalous magnetic moment > 11 x 10 “^° 

and values of /iq in the range 0.13 ± 0.05, we hnd that the point with the highest 
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tan/? 



r 

t 

h 

10 

108 (174) 

184 (306) 

132 (197) 

376 (686) 

115 (116) 

20 

154 (255) 

268 (457) 

175 (274) 

603 (990) 

116 (118) 

30 

191 (310) 

338 (560) 

212 (312) 

740 (1200) 

117 (118) 

40 

201 (340) 

357 (617) 

274 (353) 

785 (1314) 

117 (119) 

50 

208 (357) 

371 (646) 

440 (427) 

822 (1357) 

117 (119) 

55 

146 (311) 

260 (563) 

424 (676) 

606 (1237) 

115 (117) 


Table 1: Upper bounds, in GeV, on the masses of the lightest of the neutralinos, charginos, 
staus, stops and Higgs bosons for various values of tan/? if the the E821 bounds are 
imposed. The values within brackets represent the same situation when the weaker bounds 
11 X 10“^° < < 75 X 10“^° are used (see main text). 

value of mo is (in GeV) at (mo, M 1 / 2 ) = (950,300) and that with the highest value 
of Mi /2 is at (mo,Mi/ 2 ) = (600,750). The latter marks the lower end of the line 
segment of the boundary = 11 x 10 “^° which amputates the cosmologically al¬ 

lowed stripe. For the case displayed in the bottom right panel of the hgure the upper 
mass limits put on the LSP, and the lightest of the charginos, stops and the staus are 
m^ < 287,m^+ < 539,m^ < 1161,m^ < 621 (in GeV). Allowing for Ao 7 ^ 0 values, the 
upper bounds put on mo, M 1/2 increase a little and so do the aforementioned bounds on 
the sparticle masses. Thus it appears that the prospects of discovering GMSSM at a 
collider with center of mass energy ^/s = 800 GeV, such as TESLA, are not guaranteed. 
However in the allowed regions the next to the lightest neutralino, has a mass very 
close to the lightest of the charginos and hence the process e’''e“ —»■ yy', with y' subse¬ 
quently decaying to y -|- or x -|- 2 jets, is kinematically allowed for such large tan/?, 
provided the energy is increased to at least = 900 GeV. It should be noted however 
that this channel proceeds via the t-channel exchange of a selectron is suppressed due to 
the heaviness of the exchanged sfermion. 

The situation changes, however, when the strict E821 limits are imposed = 

(43.0 ± 16.0) X 10“^°. For instance in the tan/? = 55 case displayed in hgure ^ there is no 
cosmologically allowed region which obeys this bound. For the other cases, tan/? < 50, 
the maximum allowed M 1/2 is about 475 GeV, occurring at mo — 375 GeV, and the 
maximum mo is 600 GeV when M 1/2 — 300 GeV. The upper limits on the masses of the 
sparticles quoted previously reduce to m^ < 192, m^+ < 353, m/ < 775, mf < 436 all in 
GeV. However, these values refer to the limiting case Aq = 0. Scanning the parameter 
space allowing also for Aq 7 ^ 0 we obtain the upper limits displayed in the table |^. In this 
the unbracketed values correspond to the E821 limits on the 5 ^^ — 2. For completeness 
we also display, within brackets, the bounds obtained when the weaker lower bound 
> 11 X 10“^° is imposed. We see that even at TESLA with center of mass energy 
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Figure 2: In the (Mi/ 2 ,mo) plane, we display all points compatible with _ 

(43.0 ± 16.0) X 10“^° (+) and 11 x 10“^° < < 75 x 10-^° (o). All the points 

are consistent with the cosmological bound /iq = 0.13 ± 0.05 and they are grouped in 
regions, separated by dashed contours each of which is the boundary of tan/? with the 
value shown beneath. In the top region, designated by tan/? = 55, the parameter tan/? 
takes values between 50 and 55. 

a/s = 800 GeV, the prospects of discovering CMSSM are guaranteed in the ^ X~^X~ 

if the E821 bounds are imposed. 

In the figure we display in the (M 1 / 2 , fna) plane the points which are consistent both 
with the muon’s anomalous magnetic moment bounds mentioned before and cosmology, 
as well as with the other accelerators data. Each of the points is taken from a sample 
of 45,000 random points in the part of the parameter space defined by mo < 1.5 TeV, 
Mi /2 < 1-5 TeV, |Ao| < 1 TeV and 2 < tan/? < 55. All the points are consistent with 
the cosmological bound = 0.13 ± 0.05. The plus points (colored in blue) are those 
consistent with the E821 bound 27 x 10“^° < < 59 x 10“^°, while the diamonds 

(colored in green) are consistent with the more relaxed bound 11 x 10“^° < < 

75 X 10“^°. The points are grouped in regions, separated by dashed contours, each of 
which constitutes the boundary of tan/? with the value shown beneath. In the region 
designated as tan/? = 55 all points have 55 > tan/? > 50. It is seen clearly that only 
a few points in the tan/? > 50 case can survive the E821 bound. For tan/? < 50 the 
parameter M 1/2 cannot be larger than about 500 GeV, attaining its maximum value at 
mo — 400 GeV , and the maximum mo is 725 GeV occurring at M 1/2 — 275 GeV. The 
upper limits put on mo, M 1/2 result to the sparticle mass bounds displayed in the table [|. 
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Figure 3: Scatter plot of the scalar neutralino-nucleon cross section versus from a 
random sample of 45,000 points. On the top of the hgure the CDMS excluded region 
and the DAMA sensitivity region are illustrated. Pluses (+) (in blue colour) are points 
within the E821 experimental region = (43.0±16.0) x 10“^° and also cosmologically 

acceptable /iq = 0.13 ± 0.05. Diamonds (o) (in green colour) are also cosmologically 
acceptable points, but with within the region 11 x 10“^° < < 75 x 10“^°. 

Crosses (x) (in red colour) represent the rest of the random sample. The Higgs boson 
mass bound rrih > 113.5 GeV is properly taken into account. 


3 Direct Dark Matter searches 


We shall discuss now the impact of the gfj, — 2 measurements and of the Higgs mass bound 
m/j > 113.5 GeV on the direct DM searches. We are using the same random sample as 
in hgure || in order to calculate the spin-independent, X"iiucleon cross section (cr^caZar)- 
In hgure we plot the scalar x-nucleon cross section as function of the LSP mass, m^. 
On the top of the hgure the shaded region (in cyan colour) is excluded by the CDMS 


experiment . The DAMA sensitivity region (coloured in yellow) is also plotted . 


Pluses (±) (in blue colour) represent points which are both compatible with the E821 data 


a' 


SUSY _ 


^ = (43.0 ±16.0) X 10“^° and the cosmological bounds for the neutralino relic density 

^0 ~ 0-13 ± 0.05. Diamonds (o) (in green colour) are points which are cosmologically 
acceptable with respect to the aforesaid bounds, but the bound to the has been 


relaxed to its 2ct region, namely 11 < X 10^*^ < 75. The crosses (x) (in red colour) 

represent the rest of the points of our random sample. Here the Higgs boson mass bound, 
rrih > 113.5 GeV has been properly taken into account. From this hgure it is seen that 
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Figure 4: Scatter plot of the scalar neutralino-nucleon cross section versus m^, from 
a random sample of Fig. ^ Diamonds (o) (in green colour) are cosmologically accept¬ 
able points, without putting an restriction from the Crosses (x) (in red colour) 

represent points with unacceptable h^. 


the points which are compatible both the — ^ E821 and the cosmological data (crosses) 
yield cross sections of the order of 10“® — 10“® pb and the maximum value of the 
is about 200 GeV. If one allows the 2a region of the — 2 bound the lower bound of 
preferred cross sections is 10“^° pb and correspondingly the upper bound of is drifted 
up to 350 GeV. 

Gomparing hgure ^ and one can realise how g^ — 2 data constrain m-^ mass to be 
up to 200 GeV or 350 GeV for the la or 2a case respectively. In hgure H we don’t 
impose the constraints stemming from g^j^ — 2 data, therefore due to the coannihilation 
processes the cosmologically acceptable LSP mass can be heavier than 500 GeV. What is 
also important to be noticed about the direct searches of DM is that imposing the gfj, — 2 
data the lowest allowed y-nucleon cross section increased by about 1 order of magnitude, 
from 10“^^ pb to 10“^° pb. Similar results are presented in Ref. [M. This fact is very 


encouraging for the future DM direct detection experiments, with sensitivities extending 
up to 10“® pb PT[]. 


4 Conclusions 


Goncluding, we combined recent high energy physics experimental information, like the 
anomalous magnetic moment of the muon measured at E821 Brookhaven experiment and 














the light Higgs boson mass bound from LEP, with the cosmological data for DM. By 
doing so we studied the imposed constraints on the parameter space of the CMSSM and 
hence we assessed the potential of discovering SUSY, if it is based on CMSSM, at future 
colliders and DM direct searches experiments. The bounds put on the sparticle spectrum 
can guarantee that in LHC but also in a e''“e“ linear collider with center of mass energy 
^/s = 800 GeV, such as TESLA, CMSSM can be discovered. The guarantee for a linear 
collider with this energy is lost in a charged sparticle hnal state channel, if the lower 
bound on the value of 5 ^^ — 2 is lowered to its ~ 2a value, but not for the LHC. In this 
case only by increasing the center of mass energy to be ~ 1.2 TeV, a linear collider 
can End CMSSM in ff* or x^'x~ channels. The impact of the E821 experiment’s result 
along with the bound on Higgs mass is also significant for the direct DM searches. We 
found that the maximum value of the spin-independent y-nucleon cross section attained 
is of the order of 10“® pb. Moreover this cross section can not be lower than 10“^° pb, 
which is very promising for the forthcoming direct DM experiments. 
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